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Abstract. The recent discovery by Paredes et al. (2000) 
of a persistent microquasar that is positionally coincident 
with an unidentified gamma-ray source has open the pos- 
sibility that other sources in the Third EGRET Catalog 
could be interpreted as microquasars as well. In this letter 
we show that some variable unidentified EGRET sources 
in the galactic plane could be produced by faint, other- 
wise undetected microquasars with precessing jets. When 
the jet points towards the observer, gamma-ray emission 
resulting from upscattered stellar photons could be de- 
tectable yielding a variable source with weak or unde- 
| tectable counterpart at longer wavelengths. Strategies for 
£N) ■ detecting these "microblazars" with forthcoming satellites 
^ are briefly discussed. 
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1. Introduction 

The third and final EGRET catalog (Hartman ct al. 1999) 
lists about 170 gamma-ray sources not yet clearly identi- 
fied with objects detected at lower energies. Many of these 
sources concentrate towards the galactic plane and corre- 
late with the spiral arms of the Galaxy, indicating a sig- 
nificant contribution from Population I objects (Romero 
et al. 1999, Romero 2001). From a purely statistical point 
of view, at least two different populations can be iden- 
tified: one of them possibly associated with the Gould 
belt, a nearby star-forming region at ~ 600 pc, and the 
other formed by brighter sources found at lower latitudes 
(Gehrels et al. 2000, Grenier 2000). 

Several types of objects have been proposed as possi- 
ble counterparts of galactic gamma-ray sources, includ- 
ing early-type stars (both isolated and forming binary 
systems), accreting neutron stars, radio-quiet pulsars, in- 
teracting supernova remnants, and black hole candidates 
(see Romero 2001 and references therein). Very recently, 
Paredes et al. (2000), Grenier (2001) and Romero (2001) 
have suggested that microquasars can also be responsible 
for some unidentified sources. In particular, Paredes et 
al. (2000) have proposed that the microquasar LS 5039, 
a massive X-ray binary with persistent non-thermal ra- 
dio emission, is physically associated with the gamma-ray 
source 3EG J1824-1514. They suggested that the observed 
gamma-ray flux is the result of external Compton scatter- 
ing of UV photons from the high-mass stellar companion. 

The fact that there is a significant number of variable 
gamma-ray sources at low galactic latitudes (e.g. Torres 
et al. 2001, Romero et al. 2001) along with the steady 
properties of gamma-ray pulsars and interacting super- 
nova remnants, makes the idea of gamma-ray emitting mi- 
croquasars particularly attractive. Mirabel & Rodriguez 
(1999) proposed that microquasars with jets forming a 
small angle with the line of sight -by analogy with the 
unified model for AGNs- could appear as microblazars, 
namely, as sources with highly variable and enhanced non- 
thermal flux due to Doppler boosting. 
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In this letter we shall present a model for galactic vari- 
able gamma-ray sources based on the idea of precessing 
microblazars. In our model, the jet precession will be in- 
duced by the gravitational torque of the stellar companion, 
which is in a non-coplanar orbit, on the accretion disk of 
the primary (usually assumed to be a stellar black hole). 
Gamma-rays can be produced by external Compton scat- 
tering of UV stellar photons of the massive companions by 
relativistic leptons (electrons and/or positrons) far from 
the base of the jet. We shall discuss cases where the par- 
ticles arc in the form of a "blob" ejected from the central 
engine and where the particles are continuously injected 
in the form of a persistent jet. In the case of high-mass 
binaries, this mechanism for generation of gamma-rays is 
more efficient than Compton scattering upon disk pho- 
tons or self-synchrotron jet emission (Georganopoulos et 
al. 2001a) and, in combination with the gravitational ef- 
fects of the stellar companion, can render suitable vari- 
able gamma-ray sources. We shall briefly discuss how these 
sources can be identified in the light of forthcoming satel- 
lite experiments like INTEGRAL. Although microblazars 
have been previously considered in the literature (e.g. 
Aharonian & Atoyan 1998, Mirabel & Rodriguez 1999, 
Georganopoulos et al. 2001a), this is the first detailed dis- 
cussion in terms of unidentified gamma-ray sources. 



2. Variable gamma-ray emission from 
microblazars 

Let us consider a high-mass binary where accretion onto 
the compact object (a black hole) results in the produc- 
tion of twin, relativistic e + e~-pair jets propagating in op- 
posite directions. A sudden variation in the injection rate 
can result in the formation of a blob (for which we adopt 
spherical geometry) which propagates down the jet with 
a bulk Lorentz factor T. We shall follow the general treat- 
ment derived by Georganopoulos et al. (2001b) for exter- 
nal Compton scattering in cxtragalactic blazars, adapting 
it to microblazars. In the blob frame the relativistic lep- 
tons are considered to have an isotropic power-law density 
distribution n'(j ') = (fc/4-7r)7 '~ P -P(7i, 72, 7 '), where 7 ' 
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Fig. 1. Spectral high-energy distribution of the scattered 
photons for a microquasar injecting a power-law spectrum 
of electrons in the photon field of the high-mass stellar 
companion (an 07 star in this example). A cut off at 
Lorentz factors 7 ~ 10 3 has been assumed. The spectral 
energy distribution of the star is also shown (left top cor- 
ner). 



is the Lorentz factor of the leptons, k is a constant and 
■P(7i>72) 7 ') = 1 f° r 7i < 7 ' < 72, and otherwise. 

Using the Lorentz invariant n/7 2 and the relation 7 = 
D'y' between the Lorentz factor of the electrons/positrons 
in the lab frame and the blob frame respectively we get 
that the electron density in the lab frame is 
n( 7 ) = - Pp(7iA D j (1) 

47T 

where D = \T (1 — /3cos(/>)] _1 is the usual Doppler factor: 
[3 and are the bulk velocity in units of c and the viewing 
angle, respectively. 

The blob is then injected into an isotropic monoener- 
getic photon field due to the companion star (typically UV 
photons) and will then produce the inverse Compton up- 
scattering of some of these stellar photons. With respect 
to the photon field, the lab frame rate of IC interactions 
per final photon energy is, 

where <jt is the Thomson cross section. For the case of 
Thomson scattering and assuming isotropic scattering in 
the electron frame (Georganopoulos et al 2001b), 

f{x) = \{l-x)P{l/A 1 2 ,l,x), ar=j4-- (3) 
3 47^0 
The luminosity is then obtained by integrating the 
scattering rate (||) over the particle energy distribution 
(that can be deduced using the density distribution given 
in Eq. (|l|) and the volume of the blob), and multiply- 
ing by the observed photon energy em e c 2 and the num- 
ber density n p = U /eom e c 2 (where U is the energy den- 
sity of the photon field). When p is steeper than 1 and 

7m in 

to: 



L f = 



(e/4eo) 1//2 <C 7 roax = 72-D, the luminosity reduces 
kVa T c U2P- 1 



dL 
dedfl 



D 3+P 



7reo(l +?>)(3 +p) V £ o 



-(p-i)/2 



(4) 



In the Klcin-Nishima regime numerical integrations are 
required (see Georganopoulos at el. 2001b). The important 
point is that if 72 ~ 10 3 then luminosities of ~ io 36 ~ 37 erg 
s _1 can be obtained in the observer's frame at EGRET's 
energy range (100 MeV- 20 GeV). 

However, in the case of a single blob, the interaction 
time with the photon field, and hence the duration of the 
gamma-ray flare, will be strongly limited by the appar- 
ently superluminal speed of the feature. For instance, in 
the case of an 07 star and a blob with j3 — 0.98 (r = 5), 
the flux will decrease to 1/e of its original value in ~ 1 s. 
Gamma-ray production by blobs in microblazars are then 
transient phenomena difficult to detect on Earth. This is 
not the case, however, of microblazars with persistent jets. 

For a continuous relativistic flow the amplification due 
to the Doppler factor changes from D 3+p to D 2+p in 
Eq. (U) (Lind & Blandford 1985). In this case the exter- 
nal Compton scattering of stellar photons yields a stable 
source inasmuch as the flux is not perturbed. We have 
calculated the spectral energy distribution for a specific 
model with p = 2, bulk Lorentz factor L = 5, viewing an- 
gle 0=10 deg, photon field of an 07 star (photon energies 
~ 10 eV), and a high-energy cut off of 72 = 10 3 >> 71. 
The results are shown in Figure 1 , where we also show the 
energy distribution of the stellar photons. We see that lu- 
minosities of ~ 10 36 ~ 37 erg _1 can be obtained at E > 100 
MeV. Further details can be introduced through a bro- 
ken particle spectrum in orther to take into account syn- 
chrotron cooling, but these refinements are not necessary 
here. 

It is interesting to notice that whereas the gamma- 
ray emission is amplified by a factor D 2+p — _D 3+2q , the 
synchrotron jet emission results amplified only by D 2+a , 
where a = (p — 1)/2 is the synchrotron spectral index (see 
Dermer 1995). Consequently, "red" microblazars, whose 
synchrotron spectral energy distribution peaks at IR en- 
ergies, will not be ultra-luminous X-ray sources (ULXs) 
since only a small part of the radiated power goes into 
the X-ray synchrotron tail of their spectra and this emis- 
sion, in addition, is not so strongly boosted as the inverse 
Compton gamma-ray emission. Only in case of very low 
values of 72 (say around 10 2 ) we could get an X-ray source 
(Lx ~ 10 33-34 erg s _1 ) without a much more energetic 
gamma-ray source. 

The companion star in a high-mass microblazar sys- 
tem not only provides a photon field for inverse Compton 
interactions, but also a gravitational field that can exert 
a torque onto the accretion disk around the compact ob- 
ject. The effect of this torque, in a non-coplanar system, 
is to induce a Newtonian precession of the disk. If the 
jets are coupled to the disk, as it is usually thought, then 
the precession will be transmitted to them. This situation, 
which should not be confused with the geodetic precession 
(a purely General Relativity effect), has been extensively 
studied in the case of SS433 (Katz 1980) and extragalactic 
sources like 3C273 (e.g. Romero et al. 2000). 
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Fig. 2. Precessing jet model. 
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Fig. 3. Variation of the amplification factor for continuous 
jet emission as a function of time in the precessing microb- 
lazar model (for two different opening angles) . Time units 
are normalized to the precessing period. 

A sketch of the situation is presented in Figure 2. 
The disk Keplerian angular velocity is cod — (GM/rjj) 1 / 2 , 
where M is the mass of the compact object and rd is the 
radius of the precessing part of the disk. The orbital pe- 
riod is T m and the orbital radius is given by Kepler's law: 
rf n = (G/47r 2 )(m + M)T 2 , with m the mass of the star. 

The angular velocity of the tidally induced precession 



can then be approximated by ui v 
(Katz 1980, Romero et al. 2000), 



- — 3Gm cos 8 /Ar^ a uj d 
where 8 is the half- 



opening angle of the precession cone. We can now intro- 
duce a time-parametrization of the jet's viewing angle as 
(see Abraham & Romero 1999): 

(f>{t) = arcsin [x 2 + y 2 ] 1/2 (5) 

x = (cos 8 sin 4>q + sin 8 cos 4>q sin uj p t ) cos r/o — 

— sin 8 cos LOpt sin ?yo 

y = (cos 8 sin 4>q + sin 8 cos 4>q sin oj p t) sin r/o + 

+ sin 8 cos u) p t cos r] , 
where all angles are defined in Fig. 2. 

In Figure 3 we show the time evolution of the boost- 
ing amplification factor of the gamma-ray emission for the 
case of a continuous jet for two different sets of geometrical 
parameters (viewing angle of 10 deg and precession half- 
opening angle of 1 and 10 deg). Time units are normalized 
to the precession period T. We see that the flux density 
can change by a factor of ~ 6 10 3 in a single period. A 
very weak, otherwise undetected gamma-ray microblazar, 
can increase its flux due to the precession and then en- 
ter within the sensitivity of an instrument like EGRET, 
producing a variable unidentified gamma-ray source. The 
duty cycle (i.e. the fraction of time in which the source is 
highly variable) in this example is ~ 0.2T, so the source 
could appear in several EGRET viewing periods. Just to 
give a feeling of the magnitudes involved, we mention that 
for a black hole of 4 Mq, an 071a stellar companion, an or- 
bital period of 10 days, a half-opening angle of 10° for the 
precession, and an accretion disk of ~ 5 10 11 cm (ss 5 10~ 2 
a.u.), we get a precession period T ~ 100 days. The source 
could then be detectable ~ 70 days per year. This is con- 
sistent with some EGRET observations of highly variable 
sources (e.g., the case discussed by Punsly et al. 2000). 

3. Discussion 

Taking into account that there exist more than ~ 130 
high-mass X-ray binaries detected so far (e.g. Liu et al. 
2000) and that this number should be a small fraction of 
the total number of these objects in the Galaxy, it is not 
unreasonable to expect the existence of a few tens of mi- 
croblazars at mid and low galactic latitudes that could be 
responsible for the variable galactic gamma-ray sources 
detected by EGRET. Moreover, the recently discovered 
X-ray transient V4641 Sgr, which seems to harbor a ~ 9 
MQ-black hole and displays extreme superluminal veloc- 
ities, could be a microblazar (Orosz et al 2001). Since the 
stellar companion is a late B-type star, external Compton 
gamma-ray production is not expected to be very efficient 
in this case, but its high-energy emission could fall within 
GLAST sensitivity. Another interesting candidate, from 
the theoretical point of view, is the high-mass X-ray bi- 
nary LS I +61°303, which presents a one-sided jet at mil- 
liarsecond scales and evidence for a precessing accretion 
disk (Massi et al. 2001). The intrinsic jet velocity, how- 
ever, seems not to be very high: ~ 0.4c, but the source 
location is consistent with its identification with a highly 
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variable gamma-ray source 3EG J0241+6103 (Tavani et 
al. 1998, Torres et al. 2001). 

Not all microblazars should be confined to the galac- 
tic plane. Recent direct measurements of proper velocities 
in both low- (Mirabel et al. 2001) and high-mass (Ribo 
et al. 2002) microquasars show that these objects can 
present very high velocities and then could be ejected from 
the galactic plane. This is also supported by the discov- 
ery of V4641 Sgr at b ~ —4.8 deg (estimated distance: 
7.4 < d < 12.3 kpc). In the case of microquasars with 
high-mass companions, which are young objects, we could 
expect to find them up to distances ~ 100 pc (Ribo et 
al. 2002) or even more (e.g. V4641 Sgr) from the galac- 
tic plane, so they could explain some of the mid-latitude 
unidentified gamma-ray sources (Grenier 2001). A few 
low-mass microquasars could have escaped to the galac- 
tic halo and consequently might be related to some high- 
latitude unidentified EGRET detections, if the jets are 
powerful enough. 

Perhaps the best way to identified precessing gamma- 
ray microblazars is through the detection of the electron- 
positron annihilation feature in their spectra. This annihi- 
lation signature should appear as a broad, blueshifted (by 
a factor D) line in the spectrum at a few MeV, exactly 
within the energy range of IBIS imager of the forthcom- 
ing INTEGRAL satellite. Due to the precession of the jet, 
the Doppler factor will change periodically with time, and 
hence the position of the annihilation peak should oscil- 
late in energy in the lab frame around a mean value (for 
a detailed discussion of the phenomenon see Abraham et 
al. 2001). Chandra X-ray observations of non-thermal ra- 
dio sources within the EGRET location error boxes (a 
complete list of these sources is given by Torres et al. 
2001) could help to find candidates to new microquasars 
(through the detection of X-ray disk emission) , and then 
INTEGRAL exposures could be used in an attempt to 
find the annihilation line. Such a detection would be a 
remarkable discovery, since it would establish, at a same 
time, the matter content of the jets in microquasars, and 
would help to clarify the nature of some variable uniden- 
tified high-energy gamma-ray sources. 
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